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Abstract 

The Coxeter lattices, which we denote A n / m , are a family of lattices containing many of the important lattices 
in low dimensions. This includes A n , E7, E$ and their duals A* n , E$ and Eg. We consider the problem of finding 
a nearest point in a Coxeter lattice. We describe two new algorithms, one with worst case arithmetic complexity 
0(n log n) and the other with worst case complexity 0(n) where n is the dimension of the lattice. We show that 
for the particular lattices A n and A* n the algorithms reduce to simple nearest point algorithms that already exist in 
the literature. 

Index Terms 

Lattice theory, nearest point algorithm, quantization, channel coding 

I. Introduction 

The study of point lattices is of great importance in several areas of number theory, particularly the studies 
of quadratic forms, the geometry of numbers and simultaneous Diophantine approximation, and also to the prac- 
tical engineering problems of quantisation and channel coding. They are also important in studying the sphere 
packing problem and the kissing number problem [1, 2]. Lattices have recently found significant application of in 
cryptography [3, 4] and communications systems using multiple antannaes [5, 6]. 

A lattice, L, is a set of points in M™ such that 

L = {x£ K"|x = Bw,w 6 Z"} 

where B is termed the generator (or basis) matrix. We will write vectors and matrices in bold font. The ith element 
in a vector is denoted by a subscript: x t . The generator matrix for a lattice is not unique. Let M be an n x n 
matrix with integer elements such that det(M) = ±1. M is called a unimodular matrix. Then both B and BM 
are generator matrices for the lattice L. 
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Lattices are equivalent under scaling, rotation and reflection. A lattice L with generator matrix B and a lattice 
L with generator lattice B are equivalent, or isomorphic, iff 

B = aRBM 

where a > is real, R is a matrix consisting of only rotations and reflections and M is unimodular. We write 
L ~ L. 

The Voronoi region or nearest-neighbour region Vor(L) for a lattice L is the subset of M™ such that, with respect 
to a given norm, all points in Vor(i) are nearer to the origin than to any other point in L. The Voronoi region is 
an n-dimensional polytope [2]. Given some lattice point x e L we will write Vor(L) + x to denote the Voronoi 
region centered around the lattice point x. It follows that Vor(L) + x is the subset of W 1 that is nearer to x than 
any other lattice point in L. 

The nearest lattice point problem is: Given y E M™ and some lattice L whose lattice points lie in M. n , find a lattice 
point x e L such that the Euclidean distance between y and x is minimised. We use the notation NearestPt(y, L) 
to denote the nearest point to y in the lattice L. It follows from the definition of the Voronoi region that 1 

x = NearestPt(y, L) <^ y e Vor(L) + x 

The nearest lattice point problem has significant practical application. If the lattice is used for vector quantisation 
then the nearest lattice point corresponds to the minimum-distortion point. If the lattice is used as a code for 
a Gaussian channel, then the nearest lattice point corresponds to maximum likelihood decoding [7]. The closely 
related shortest vector problem has been used in public key cryptography [3, 4, 8, 9, 10]. 

Van Emde Boas [11] and Ajtai [12] have shown that the nearest lattice point problem is NP-complete under 
certain conditions when the lattice itself, or rather a basis thereof, is considered as an additional input parameter. 
It has even been shown that finding approximately nearest points is NP-hard [8, 13, 14]. Nevertheless, algorithms 
exist that can compute the nearest lattice point in reasonable time if the dimension is small [15, 16, 17]. One such 
algorithm introduced by Pohst [17] in 1981 was popularised in signal processing and communications fields by 
Viterbo and Boutros [16] and has since been called the sphere decoder. 

For specific lattices, the nearest point problem is considerably easier and for many classical lattices, fast nearest 
point algorithms are known [1, 2, 7, 18, 19, 20]. 

The Coxeter lattices, denoted A n / m , are a family of lattices first described by H.S.M. Coxeter [21, 22]. 



4 /m = {Qx|xeZ n+1 ,x'lmodra = 0} (1) 
where Q is the orthogonal projection matrix 

«-('-;£)■ 

'There is a slight technical deficiency here. We actually require to define half of the faces of Vor(L) to be closed and half to be open. Ties 
in NearestPt(y, L) can then be broken accordingly. 
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I is the (n + 1) x (n + 1) identity matrix, 1 = [1, 1, 1, . . .]' and ' indicates the vector or matrix transpose. If m 
does not divide n + 1 then A n / m — A n j X . Hence, in the sequel, we assume that m divides n + 1. 

A simple geometric description of A n / m is to consider the subset consisting of the points of Z n+1 whose 
coordinate-sum is divisible by to. This subset consists of points that lie in 'layers' parallel to the hyperplane 
orthogonal to 1. By projecting the subset orthogonally to 1 we obtain a set of points equivalent to the n-dimensinal 
lattice A n/m . 

The family of Coxeter lattices contains many of the important lattices in low dimension. The family is related 
to the well studied root lattice A n and its dual lattice A* n . When m = 1 

A n/1 = A* n = {Qx | x e Z" +1 } (3) 

and when to = n + 1 

A n/n+ i = A n = {x e | x'l = 0} (4) 
It follows that A n C A n / m C A* n [2, 22]. Note that A n / m C A„/ fc whenever k < m and therefore 

Vor(A„ /fe ) C Vor(A„ /m ). (5) 

Other isomorphisms exist: A s / 3 ~ _E 8 ~ £?g, A 7 / 4 ~ £7 and A 7 / 2 — £7. Of significant practical interest is 
the lattice E s ~ ^8/3- Due to its excellent packing and quatising properties E s has found applications to trellis 
codes [23, 24, 25, 26] and vector quantisation [2, 27, 28]. The particular representation of E s as A s / 3 was used 
by Secord and deBuda to create a code with a spectral null at DC [29]. 

The lattice A* n ~ is also of practical interest. It gives the thinnest sphere-covering in all dimensions up to 
8 [2] and has found application in a number of estimation problems including period estimation from sparse timing 
data [30, 31], frequency estimation [32], direction of arrival estimation [33] and noncoherent detection [34]. 

The paper is organised as follows. Section II describes a log -linear- time nearest point algorithm for A n / m . This 
algorithm is a generalisation of a nearest point algorithm for A* n that was derived in [19]. Section III improves this 
to worst case linear-time. The speedup employs both a partial sorting procedure called a bucket sort [35] and also 
the linear-time Rivest-Tarjan selection algorithm [36, 37, 38, 39]. In Section IV we show how the discussed nearest 
point algorithms for the Coxeter lattices reduce to simple nearest point algorithms for A n and A* n that already exist 
in the literature [2, 19, 20]. In Section V we review a simple nearest point algorithm for A n / m based on translates 
of the lattice A n . This algorithm was previously described by Conway and Sloane [7, 18] but not directly applied 
to the Coxeter lattices. The algorithm requires 0(n 2 ) arithmetic operations in the worst case. In Section VI we 
evaluate the practical computational performance of the algorithms. 

II. Log-linear-time algorithm 

In this section we describe a nearest point algorithm for A n / m that requires 0(n log n) operations in the worst 
case. This algorithm is a generalisation of the nearest point algorithm for A* n described in [19]. To describe the 
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algorithm we first require to derive some properties of the Voronoi region of A n / m . This is done in Lemmata 1 
and 2. We firstly require the follow definitions. 

Let H be the hyperplane in K" +1 orthogonal to 1. H is typically refered to as the zero-mean-plane. For some 
lattice L we will use the notation Vor#(L) to denote the region Vor(L) n H. For example Vor#(A„) is the 
crossection of Vor(A„) lying in the hyperplane H. Given some region R C H we define the n-volume of R as 
voLy(-R). For example, the n-volume of Vor#(A„) is denoted by vol#(Vor# (A„)). 

Given a set of n-dimensional vectors S and suitable matrix M we will write MS* to denote the set with elements 
Ms for all s £ S. For example Q Vor(Z™ +1 ) denotes the region of space that results from projecting Vor(Z™ +1 ) 
onto the hyperplane H. 

Lemma 1. 

QVor(Z" +1 ) C Yor H (A n ) 

Proof: Let y G Vor(Z" +1 ). Decompose y into orthogonal components so that y = Qy + tl for some tel. 
Then Qy e Q Vor(Z" +1 ). Assume that Qy f Vor H (A n ). Then there exists some x E A n such that 

||x-Qy|| 2 < ||0-Qy|| 2 ||x-y + tl|| 2 < ||y-tl|| 2 

^||x-y|| 2 + 2tx'l<!|y|| 2 . 

By definition (4) x'l = and so ||x-y|| 2 < ||y|| 2 . This violates that y e Vor(Z™ +1 ) and hence Qy e Yor H (A n ). 

■ 

Lemma 2. 

Vor ff (A n/m )CQVor(Z n+1 ) 

with equality only when m = n + 1. 

Proof: When m = n + 1, A n / n+ i = A n . The n-volume voltf(Vor#(^4„)) = \/n + 1 [2]. From Berger et al. 
[40] we find that the n-volume of the projected polytope vol#(Q Vor(Z™ +1 )) = \/n + 1 also. As Vor ff (^4 n ) and 
Q Vor(Z rl+1 ) are convex polytopes it follows from Lemma 1 that 

Vor H (A n ) = QVor(Z" +1 ). 

The proof follows from the fact that Vorjj(^4 n / TO ) C Vor h(A u ) for all m (5). ■ 
We will now prove Lemma 3 from which our algorithm is derived. We firstly need the following definition. 
Given two sets A and B we let A + B be their Minkowski sum. That is, x e A + B iff x = a + b where a E A 
and b e B. We will also write IK to denote the line of points lr for all rel. Then Vor ii{A n / m ) + 1M is an 
infinite cylinder with cross-section Vor#(A„/ m ). It follows that VovH(A n / m ) + 1R = Yov{A n / m ) 

Lemma 3.1fx = Qk is a closest point in A n j m to y S R™ +1 then there exists some A e R for which k is a 
closest point in Z" +1 to y + Al. 
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Proof: As Qk is the nearest point to y then for all A e R 

y + 1A G Voi(A n/m ) + Qk = Vor H {A n/m ) + k + 1R. 
It follows from Lemma 2 that 

Vor„(A n/m ) + k + 1R C Q Vor(Z Il+1 ) + k + 1R. 
Then y + 1A e Q Vor(Z" +1 ) + k + 1R and for some A e R 

y + lAe Vor(Z" +1 )+k 

The proof now follows from the definition of the Voronoi region. ■ 
Now consider the function f :ln Z n+1 defined so that 

f (A) = Ly + Al] (6) 

where [•] applied to a vector denotes the vector in which each element is rounded to a nearest integer 2 . That is, 
f (A) gives a nearest point in Z™ +1 to y + Al as a function of A. Observe that f (A + 1) = f (A) + 1. Hence, 

Qf(A + 1) = Qf(A). (7) 

Lemma 3 implies there exists some A e R such that x = Qf(A) is a closest point to y. Furthermore, we see 
from (7) that A can be found within an interval of length 1 . Hence, if we define the set 

5 = {f(A) | A e [0,1)} 

then QS contains a closest point in A n / m to y. In order to evaluate the elements in S we require the following 
function. 

Definition 1. (sort indices) 

We define the function 

s = sortindices(z) 

to take a vector z of length n + 1 and return a vector s of indices such that 

%S\ ^ z S2 i> z S3 ^ • • • > z Sn+1 

Let 

s = sortindices({y}) 

2 The direction of rounding for half-integers is not important so long as it's consistent. The authors have chosen to round up half-integers in 
their own implementation. 
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where {g} = g — [g] denotes the centered fractional part of g e M and we define {•} to operate on vectors by 
taking the centered fractional part of each element in the vector. It is clear that S contains at most n + 2 vectors, 
i.e., 

{[y],[y] +e ai ,Lyl + e Sl +e S2 ,..., 

Ly] +e Sl +--- + e s „ +1 } (8) 

where ei is a vector of O's with a 1 in the ith position. It can be seen that the last vector listed in the set is simply 
Ly] + 1 and so, once multiplied by Q, the first and the last vector are identical. 
We can define the set W C S such that 

W = {x e S | x • 1 mod to = 0}. (9) 

Noting (1) then QW contains the nearest point in A n / m to y. 

An algorithm suggests itself: test each of the distinct vectors in QW and find the closest one to y. This is 
the principle of the algorithm we propose in this Section. It remains to show that this can be done in 0(n log n) 
arithmetic operations. 

We label the elements of S according to the order given in (8). That is, we set uo = and, for i = 1, . . . , n, 

Ui = Ui_i+e Si . (10) 

Let Zi = y — Ui. Clearly, Zo = {y}. Decompose y into orthogonal components Qy and tl for some fgl. The 
squared distance between Qu^ and y is 

||y-Qui|| 2 = di + t 2 (n+l) (11) 

where we define di as 

d l = HQzjII 2 



z I 

-1 



ag. (i2, 

n + 1 



n+l 

We know that the nearest point to y is that Qu^ such that Uj e W which minimizes (11). Since the term t 2 (n + 1) 
is independent of the index i, we can ignore it. That is, it is sufficient to minimize di, i = 0, . . . , n. 

We now show that di can be calculated inexpensively in a recursive fashion. We define two new quantities, 
di = z-1 and /3j = z-Zj. Clearly di — fli — a 1/n+\. From (10), 

oti = z-1 = (zj_i - e Si )'l = £*,_!- 1 (13) 

and 

Pi = z'^i = (zj_i - e Si )'(zj_i - e Si ) = Pi-i - 2{y Si } + 1. (14) 

Algorithm 1 now follows. The main loop beginning at line 8 calculates the on and fa recursively. There is no 
need to retain their previous values, so the subscripts are dropped. The variable D maintains the minimum value of 
the (implicitly calculated values of) di so far encountered, and k the corresponding index. The variable 7 maintains 
the value of u-1 mod m which must equal in order for <G W. 
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Each line of the main loop requires O(l) arithmetic computations so the loop (and that on line 15) requires 
0(n) in total. The function sortindices(z) requires sorting n+1 elements. This requires 0(n log n) arithmetic 
operations. The vector operations on lines 2-5 all require 0(n) operations and the matrix multiplication on line 17 
can be performed in 0(n) operations as 

l'u 

Qu = u 1. 

n+1 

It can be seen, then, that the computational cost of the algorithm is dominated by the sortindices(-) function and 
is therefore 0(n log n). 

This algorithm is similar to the nearest point algorithm for A* n described in [19]. The significant difference is the 
addition of 7 = on line 9. This ensures that the lattice points considered are elements of A n / m i.e. they satisfy 
(1). We further discuss the relationship between the algorithms in Section IV. 

Input: y e W l+1 

1 u = Ly] 

2 z = y — u 

3 a = z'l 

4 (3 = z'z 

57 = u'l mod m 

6 s = sortindices(z) 

7 D = 00 

8 for i = 1 to n + 1 do 



9 

10 
11 

12 

13 
14 



if /3 - < D and 7 = then 

k = i-l 

a = a — 1 

f3 = f3- 2z Si + 1 

7 = (7 + 1) mod m 



15 for i = 1 to k do 



16 



u Si = u s . + 1 



17 x = Qu 
is return x 

Algorithm 1: Algorithm to find a nearest lattice point in A n / m to y € R n+1 that requires 0(n log n) arithmetic operations 



III. Linear-time algorithm 

In the previous Section we showed that the nearest point to y in A n / m lies in the set QW (9). We will show that 
some of the elements of QW can be immediately excluded from consideration. This property leads to a nearest 
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point algorithm that requires at most 0(n) arithmetic operations. 
Lemma 4. Suppose, for some integers i, m > 0, k > 2, that 

{Vsi} ~ {Vs i + km \ < n J r \' 

Then the minimum of the d i+cm , c = 0, . . . , k, occurs at c = or c = k. 

Proof: The proof proceeds by contradiction. Suppose, to the contrary, that 

di-\-cm ^ di and di-\- C m < ^H^m • 

Observe that 

2 cm 



2a,- cm — (cm) v — • 
d i+cm ~di= + - 2 ^+, »• 



Now, since {j/ s<+i } < {j/s 4 }> it follows that 



2cticm — (cm) 



2 



di+cm - di > n^Ti h cm ^ _ 2 {y^}) 



n + 

With the assumption that d i+cm — di < 0, we have that 

2a; — cm 

— < 2{y s .} - 1. 

n+1 iy s/ 

Similarly, observe that 

_ 2cti(k - c)m - (A: 2 - c 2 )m 2 
n+1 

km 

+ (i-2{y Sl+ J). 

j=cm+l 

Since > {y Si+km }, it follows that 

2ai(k — c)m — (k 2 - c 2 )m 2 



n+1 

+ (k-c)m(l-2{y Si+km }). 



With the assumption that di + k m ~ d i+cm > 0, we have that 

2oti — cm km „ , , 

Equations (16) and (17) together imply that 

km 



{y Si } - {y St+km } > 



2(n + l)' 
which contradicts (15) because k > 2. 

From S we can construct the following q = ™+ 1 /m subsets 

'm(j - 1) mj 



Uj = { u 2 | 0.5 - {y s J e 



n+1 n+1 
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where j — 1, • • • ,q. Note that QS = Q Uj=i ^i- We 316 interested in the elements of UjtlW. Let g be the smallest 
integer such that u g e Uj n W- 7 . Let p be the largest integer such that u p e [/,■ n W. It follows that p = g + km 
for some fc g Z. Also, from (18) 

m 

{l/» 9 }-{lfep} < ^-py- 

It then follows from Lemma 4 that (11) is minimised either by u s or u p and not by any e Uj D W where 
g < i < p. We see that for each set QUj there are at most two elements that are candidates for the nearest point. 
An algorithm can be constructed as follows: test the (at most two) candidates in each set QUj and return the closest 
one to y. We will now show how this can be achieved in linear time. 
We construct q sets 



Bj = \i | 0.5 - { Vi } e 



m(j - 1) mj 



n+l n+1 



(19) 



and the related sets 



It follows that 



Definition 2. (quick partition) 

We define the function 



Kj = (jB t . 
t=i 



teK, 



b = quickpartition(z, Bj, c) 



fo take a vector z a«<i integer c= 1, . . . , \Bj\ and return a vector b of length \Bj\ such that for i = 1, . . . ,c — 1 
one/ t = c + 1, . . . , |.Bj| 

Somewhat surprisingly quickpartition(z, Bj, c) can be implemented such that the required number of operations 
is 0(\Bj\). This is facilitated by the Rivest-Tarjan selection algorithm [36, 37, 38, 39]. We can compute 

b = quickpartition(z, Bj , c) (20) 

for some integer 1 < c < | Bj \ . Then 

c 

U|^_ 1 | +c = U| X ._ 1 |+^e 6t . (21) 
tei 

Let g be the smallest integer such that 1 < g < \Bj\ and 

1 • ui^.^+g mod m = (22) 
and let p be the largest integer such that 1 < p < \Bj\ and 

1-U|^_ 1 | +P modm = 0. (23) 
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From the previous discussion the only candidates for the nearest point out of the elements 

Q {uik^i+1, . . . ,U| Kj l | + | B3 |} = QUj 

are Qu|#- | +g and Qu|jf | +p . We can compute these quickly using the quickpartition(-) function as in (20) 
and (21). 

Algorithm 2 now follows. Lines 2-5 construct the sets Bj. The main loop on line 12 then computes the values 
of g and p for each Bj. We define the function 

b = quickpartition2(z, Bj,g,p) 

to return b so that for i = 1, . . . , g — 1 and t = g + 1, . . . ,p — 1 and c = p + 1, . . . , \Bj \ 

Zbi > Zb g > Zb t > z bp > Zb c . 

Notice that quickpartition2(-) can be performed by two consecutive iterations of the Rivest-Tarjan algorithm and 
therefore requires 0(\Bj\) operations. The d\K } \+ g and d\ K .\ +p are computed within the loop on line 16 and the 
index of the nearest lattice point is stored using the variable k* . The concatenate(w, b) function on line 24 adds 
the elements of b to the end of the array w. This can be performed in 0(\Bj\) operations. Lines 25-27 recovers 
the nearest lattice point using w and k*. 

In practice the Bj can be implemented as a list so that the set insertion operation on line 5 can be performed in 
constant time. Then the loops on lines 2 and 3 require 0(n) arithmetic operations. The operations inside the main 
loop on line 12 require 0(\Bj\) operations. The complexity of these loops is then 

n+l/ m 

0(\Bj\) = 0(n) 

3 = 1 

The remaining lines require 0{n) or less operations. The algorithm then requires 0(n) arithmetic operations. 

IV. Specific algorithms for A n and A* n 

For the lattices A n = A n / n+l and A* n = A n ^ Algorithms 1 and 2 reduce to simpler algorithms that have 
previously been described in the literature. For A n a log-linear time algorithm similar to that of Conway and 
Sloane [7, 18] is derived from Algorithm 1 by noting that only one iteration in the main loop on line 8 will satisfy 
7 = 0. Algorithm 3 now follows. 

A simple linear-time algorithm for A n can be constructed from Algorithm 3 by replacing the sortindices(-) 
function on line 2 with quickpartition(-). Pseudocode is provided in Algorithm 4. In effect this is a modification 
of Algorithm 2 where the sets from (19) are replaced by the single set {1, 2, . . . , n + 1}. This algorithm has 
previously been suggested by A. M. Odlyzko [2, page 448]. 

For A* n a log-linear time algorithm identical to that described in [19] can be derived from Algorithm 1 by 
noting that 7 mod 1 = for all 7. A linear-time algorithm for A* n can be constructed from Algorithm 2 by noting 
that g — 1 (22) and p = \Bj\ (23) for all Bj where j = l,2,...,n+l. This removes the need for using the 
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Input: y e 

1 z = y - |_y] 

2 for j = 1 to q do Bj 

3 for i = 1 to n + 1 do 



j = q- [q( Zl + l/ 2 )J 



e u = iyi 

7 Q = z'l 

8 /3 = z'z 

97 = u'l mod m 

10 fc = 1 

11 _D = 00 

12 for j = 1 to g do 



13 

14 
15 
16 
17 
18 
19 
20 
21 
22 

23 

24 



<; = m — 7 

p = l-Bjl — (l-Bjl +7) mod m 
b = quickpartition2(z, Bj, g,p) 
for i = 1 to I -Bj j do 

a = a — 1 
/? = /?- 2z bt + 1 
7 = (7 + 1) mod m 

if (i — g or i = p) and /3 — a2 /n+i < D then 
D = 0- a 2 /(n + 1) 
fc* = fc 

fc = fc + l 

concatenate(w, b) 



25 for i = 1 to fc* do 



26 



27 X = Qu 

28 return x 

Algorithm 2: Algorithm to find a nearest lattice point in A n / m to y G K n+1 that requires 0(n) arithmetic operations 
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Input: y E R n+1 

1 7 = (n + 1 — |_y~l'l) m °d n + 1 

2 s = sortindices({y}) 

s u = [y] 

4 foreach i = 1 to 7 do 

5 u Si = u Si + 1 

6 x = Qu 

7 return x 

Algorithm 3: Algorithm to find a nearest lattice point in A„ to y € R™ that requires 0(n log n) operations 

Input: y e M. n+1 

1 7 = (n + 1 — |_y| 1) mod n + 1 

2 b = quickpartition({y} , {1, 2, . . . , n + 1}, 7) 

3 u = |_y] 

4 for i = 1 to 7 do 

5 u b% = u hi + 1 

6 x = Qu 

7 return x 

Algorithm 4: Algorithm to find a nearest lattice point in A n to y € K n that requires O(n) operations 

quickpartition2(-) function. A further simplification is noted in [20] where it was shown that the nearest point is 
one of the Qii|x-| where j = 0, • • • , n. The reader is referred to [20] for further details. The proofs used in [20] 
are significantly different to those in this paper and are only applicable to A* n . 

V. Algorithm based on glue vectors 

In this section we describe a simple nearest point algorithm for A n / m . This algorithm was described by Conway 
and Sloane [7, 18] but not directly applied to the Coxeter lattices. The algorithm has worst case complexity 0{n 2 ). 
A n / m can be constructed by gluing translates of the lattice A n [2]. That is 

9-1 

A n /m = IJ ([im] +A n ) (24) 

i=0 

where q = ™+ 1 /m and [i] are called glue vectors and are denned as 

[i] = ^g-j- {i, ■ ■■ , ■ ■■ , -j) (25) 

j times i times 

for i G {0, . . . , n} with i + j = n + 1. Following the notation of Conway and Sloane the glue vectors will not be 
written in boldface. Instead they are indicated by square brackets. 
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Noting that A n / m can be constructed as a union of q translates of the lattice A n we can use a nearest point 
algorithm for A n to find the nearest point in each of the translates. The translate containing the closest point yields the 
nearest point in A n / m . A pseudocode implementation is provided in Algorithm 5. The function NearestPt(y, A n ) 
can be implemented by either Algorithm 3 or 4 of Section IV. 

Input: y e M. n 

1 f = 00 

2 for i = to q 1 do 

3 x = NearestPt(y — [im], A n ) + [im] 

4 if ||x-y|| < D then 

5 Xnp = X 

L^ = ||x-y|| 
7 return x NP 

Algorithm 5: Nearest point algorithm for A n/m using glue vectors 

The algorithm requires iterating NearestPt(y, A n ) q times. Assuming that NearestPt(y, A n ) is implemented 
using the linear time algorithm (Algorithm 4) then if q is a constant this yields a linear-time algorithm. At worst 
q may grow linearly with n. In this case the algorithm requires 0(n 2 ) operations. 

VI. Run-time analysis 

In this section we tabulate some practical computation times attained with the nearest point algorithms described 
in Sections V, II and III and also some of the specialised algorithms for A n and A* n discussed in Section IV. The 
algorithms were written in Java and the computer used is a 900 MHz Intel Celeron M. 

Table VI shows the computation times for the three algorithms from Sections V, II and III for the lattice A n / 4 and 
q = n + 1 /4. It is evident that the linear-time algorithm is the fastest. The glue vector algorithm is significantly slower 
for large n. By comparison, Table VI shows the computation times for the algorithms with A n / m for m = n + 1 /4 and 
q = 4. The glue vector algorithm now performs similarly to the other algorithms. This behaviour is expected. As 
discussed in Section V the glue vector algorithm has linear complexity when q is constant, but quadratic complexity 
when q increases with n. 

Tables VI and VI show the performance of the linear-time Coxeter lattice algorithm compared to the specialised 
algorithms for the lattices A* n and A n discussed in Section IV. It is evident that the specialised algorithms are faster. 
This behaviour is expected as the specialised algorithms have less computational overhead. 

VII. Conclusion 

In this paper we have described two new nearest point algorithms for the Coxeter lattices. The first algorithm is a 
generalisation of the nearest point algorithm for A* n described in [19] and requires 0(n log n) arithmetic operations. 
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TABLE I 

Computation time in seconds for A n / 4 for 10 5 trials 



Algorithm 


n=25 


n=100 


n=1000 


O(n) 


6.14 


18.89 


165.57 


0(n log n) 


6.83 


21.51 


205.36 


0(n 2 ) 


13.66 


161.80 


> 10 4 



TABLE II 

COMPUTATION TIME IN SECONDS FOR A n / m WITH m = 2±i FOR 10 5 TRIALS 

Algorithm n=25 n=100 n=1000 

O(n) 6.67 17.78 157.66 

0(n log n) 21.33 9.27 209.23 
0(n 2 ) 10.71 35.24 317.14 



The second algorithm requires 0(n) operations in the worst case. The second algorithm makes use of a partial 
sorting procedure called a bucket sort [35] and also the linear-time Rivest-Tarjan selection algorithm [36, 37, 38, 39]. 
In Section IV we showed how the log-linear and linear-time algorithms for the Coxeter lattices reduce to simple 
nearest point algorithms for A n and A* n that already exist in the literature [2, 19, 20]. 
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